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ABSTRACT 

The main goal of this work is to better understand foam behavior both on the Earth and in 
microgravity conditions and to determine the relation between a foam’s structure and wetness 
and its rheological properties. Our experiments focused on the effects of the bubble size 
distribution (BSD) on the foam behavior under gradual or stepwise in the liquid flow rate and on 
the onset of the convective instability. We were able to show experimentally, that the BSD 
affects foam rheology very strongly so any theory must take foam texture into account. 
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Figure 1: Behavior of stabilized foams as a function of liquid flow rate q. The BSD histograms 
corresponding to data sets 1-3, 4-6 and 7-10 are shown in a lower row. 

When we change liquid flow rate by small increments, we detected four different stages 
of foam evolution. At low flow rates, increase in liquid supply produces only thickening of 
membranes and borders but does not induces bubbles motion (Stagel). For larger flow rate, but 
before the onset of global convection, we observed two possible types of foam behavior. One of 
them includes localized bubbles rearrangements in response to variations in liquid supply, which 
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stop after few minutes (Stage 2). For almost monodispersed foams (data sets 4-6) this stage 
presents for the narrow range of liquid flow rates or may be even absent. If the BSD is wide, 
with appreciable fraction of large bubbles (data sets 7-10), these rearrangements are typical for 
larger liquid flow rates range. Such “placing right bubbles to the right place” reduces stress in 
foam lattice and delays the onset of global convection for polydispersed foams. This stage was 
not observed for foams containing relatively large fraction of small bubbles. Instead, we have 
seen “scattered motion” (Stage 3): bubble cluster start to move, stops after few seconds or 
minutes, then similar motion begins at other place at the foam column, and so on. Such flow 
pattern persists as long as liquid flow rate does not change. Both Stages 2 and 3 for the same 
foam were never observed. “Scattered” motion always starts from the displacement of the 
smallest bubbles while the larger ones become involved later. Small bubbles, which may travel 
along the borders between large ones, may acts as increase effective viscosity of basic fluid, such 
that the passage of small bubbles between larger ones locally destabilizes the foam. The large- 
bubble motion stops after structure adjustment, but may start again if another cluster of small 
bubbles will be provided by incoming liquid flow. 

We also examine the development of convective pattern under rapid increase in liquid 
flow rate from subcritical (Stage 1) to supercritical value (Stage 4). Typically, the melting wave 
propagates down along one side of the cell while the much slower upward foam flow develops 
on the other side. During few seconds after the flow rate change, there is no bubble motion while 
liquid flow through the foam slightly increases. After that, the entire foam downstream of the 
convective front moves down. Bubbles slightly compresses, but still do not moves relative to 
each other. When the melting wave passes few centimeters, the upward counterflow begins to 
develop. If the difference in flow rate is large, bubbles at the front of the wave may be 
significantly compressed. Behind the front, the borders thicknesses increase and bubbles may 
rotate while moving down. When the melting front reaches the bottom of the column, steady 
convective pattern develops. We captured the foam motion by digital video camera and measured 
the averaged image density. Borders are darker than lamellae, such that the average density 
provides information about liquid distribution in the foam. During the initial prewetting, the 
wetness of motionless foam increases uniformly over entire column. For the compression stage, 
we see the growth of wetness fluctuations and beginning of development of wetness maximum in 
front of approaching melting wave. The position of moving melting front is clearly seen on the 
vertical density profile. The wetness difference behind the front and downstream and the front 
slope are related to the wave speed and, in turn, to jump in liquid flow rates. When the 
convective roll is fully developed, vertical profiles of relative density again become uniform, 
while on horizontal profiles typically appear two maxima corresponding to downward and 
upward foam flows. The amplitude of these intensity maxima is roughly proportional to the foam 
velocity. The surprising is the overall decrease of the relative density during the compression 
stage. The physical mechanism behind it is not clear and requires additional study. 
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Convection in “Stabilized” Foams 



standing foam 


convection starts 
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Stepwise Increase in Liquid Supply 


• Compression wave propagates downward, preceding 
melting wave. 

• Fast increase in thickness of Plateau membranes and 
Plateau borders. 

• Second increase in thickness and adjustment of 
curvatures of the membranes and Plateau borders. 

• The degree of compression is roughly proportional to the 
increase in liquid flow rate. 

• Convective roll develops. 
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Propagation of the Melting Front. 
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The liquid flow rate increased from 0.375 to 2.03 cm 3 /sec. (Right) 
Vertical and horizontal profiles of relative image density. * - initial 
prewetting, no motion; 0 - compression, O - melting front passes. 
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Critical Flow Rate Measurements 

• We increased the flow rate slowly in small increments until the 
foam began to move. We defined the critical flow rate as the flow 
rate at which global convection, which affects every bubble, starts. 

• We found four stages of foam evolution under the gradual increase 
in liquid supply: 

1. Membrane thickness increases, no bubble motion. 

2. Localized bubble rearrangements. The motion stops when “right” 
bubbles are placed into “right” places. 

3. Scattered movement in different areas of foam column. 

4. Global convection. 

• Critical flow is measured as a function of viscosity of basic fluid, 
BSD and fluid fraction for foams in cylindrical and rectangular (6.2 
cm x 4.7 cm x 25 cm) foam cells. 
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Effect of Bubble Size Distribution - 1 
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Effect of Bubble Size Distribution - 2 


rectangular cell 



average bubble diameter d, mm 
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Effect of Bubble Size Distribution - 3 


• Foams with smaller average bubble diameter are more 
UNSTABLE than those with larger average bubble 
diameter. 

• Homogeneous foams are more UNSTABLE to global 
convection than polydisperse foams with the same average 
bubble diameter. Additional stress may appear within the 
foam lattice because uniform bubbles are not able to fill 
the cell perfectly. 

• Foams are very STABLE when average bubble diameter is 
comparable with cell dimensions. 

• The results for small bubbles are hardly reproducible. 
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Effect of Base Fluid Viscosity 

cylindrical cell 



• Data scattering is mostly due to variations in bubble size 
distribution (BSD). Critical flow rate becomes almost independent 
of viscosity at p > 10 cp 

• Only one convective roll is found in our cylindrical foam cell 
(diameter ~ 4.5 cm) if the flow rate is high enough. 
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Critical Liquid Fraction 


• In the rectangular 
foam cell (6.2 cm x 
4.7 cm), one or two 
rolls may be formed. 

• Presence of cell 
corners reduces foam 
stability and increases 
its sensitivity to 
nonuniformities in 
liquid supply. 

'0 1 2 3 4 5 6 7 8 9 10 11 

average bubble diameter d, mm 

• In stressed foams, the critical wetness is significantly below the 
close packing loss transition (0 C ~ 0.36). 
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The estimate for average 
liquid film thickness h 


Scaling 

For one draining soap film with parabolic 
velocity profile inside: 


For one bubble within foam lattice: 
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Global convection starts when border 
curvatures cannot be adjusted to pass 
large incoming liquid flow. 
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Conclusions 


• Stressed foams can undergo melting transition when foam 
wetness is significantly below the one which corresponds to 
the close packing loss transition. 

• This melting transition may start when the balance between 
border curvatures (local pressure gradient) and local flow in 
membranes and boarders breaks down. 

• The critical parameters for this transition depend heavily not 
only on the bubble size distribution but also on the spatial 
distributions of bubbles. 
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